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Numerical  Simulation  of  the  Electromagnetic  Instability  of  an 
Intense  Beam  in  a  Quadrupole  Focusing  System 


I.  Introduction 

Discrete  quadrupole  focusing  systems,  also  called  FODO  lattices  (for  focusing  field, 
zero  field,  defocusing  field,  zero  field),  have  been  used  to  transport  charged  particle  beams 
for  a  variety  of  applications.  A  FODO  lattice  is  an  alternative  to  helical  quadrupole 
windings  for  strong  focusing.  Helical  quadrupole  (stellarator)  focusing  systems  have  been 
found  to  be  subject  to  an  electromagnetic  instability,  which  we  referred  to  as  the  three- wave 
instability.^-®  In  the  helical  quadrupole  case,  the  instability  can  be  extremely  disruptive  to 
beam  transport.  This  has  been  observed  in  numerical  simulations^’®  and  in  experiments.^ 

Here,  we  consider  the  corresponding  instability  for  transverse  perturbations  of  an  elec¬ 
tron  beam  interacting  with  a  FODO  lattice  field  and  a  TEi\  waveguide  mode.  Specifically, 
we  consider  an  alternating  gradient  quadrupole  field  (5,*,^,,,),  where 

B^j.  =  -BqkqCOs{nkgz)y,  Bgy  =  ~BqkgCOs{nkqz)x,  (la  -  b) 

Bgkg  is  the  peak  field  gradient,  kg  =  27r/Aq  and  A,  is  the  period  of  the  quadrupole  field. 
This  configuration  has  been  studied  theoretically®’*  and,  unlike  the  helical  quadrupole 
case,  a  useful  stable  regime  does  not  appear  to  exist.®  Our  present  numerical  study  was 
motivated  by  the  apparent  need  to  discover  a  detuning  mechanism  for  this  instability, 
perhaps  analogous  to  that  found  in  Ref.  3  for  the  heliceil  quadrupole  configuration. 

In  Ref.  5,  the  dispersion  relation  for  the  instability  was  derived  and  numerical  solutions 
were  presented.  These  were  shown  to  be  in  good  agreement  with  approximate  analytical 
expressions  for  the  spatial  and  temporal  growth  rates.  It  was  also  shown  that  for  given 
parameters,  peak  growth  rates  occur  at  a  single  frequency  but  with  multiple  wave  numbers. 
Futhermore,  it  was  found  that  secondarv  modes  also  exist,  with  multiple  wave  numbers  at 
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each  frequency.  Numerical  solutions  to  the  dispersion  relation  show  that  these  secondary 
modes  have  significantly  lower  growth  rates  than  the  primary  mode.  Finally,  it  was  shown 
that  the  primary  mode  is  “cut-off”  when 

Y  +  -  -  2kgKg  <  y  j ,  (2) 

where  Kq  =  =  eBq/'ym^c,  vi,  is  the  beam  velocity,  7  =  (1  —  v\l<?)^l^ ,e  is  the 

elementary  charge  (assumed  positive),  is  the  electron  mass,  c  is  the  speed  of  light,  finTg 
is  the  first  postive  zero  of  the  Bessel  function  Tg  is  the  radius  of  the  waveguide  through 
which  the  beam  is  being  transported,  and  ’  indicates  a  derivative  of  J\  with  respect  to  its 
argument. 

As  an  example,  an  (a;,  k)  diagram  is  given  in  Fig.  1  for  a  case  in  which  the  quadrupole 
gradient  is  Bqkq  =  1000  G/cm,  the  wave  number  is  kq  =  0.55  cm~^,  the  waveguide  radius 
Tg  =  3  cm  and  7  =  7.  For  clarity,  the  diagram  corresponds  to  the  limit  of  zero  beam 
current.  In  Fig.  1,  the  primary  mode  occurs  at  points  indicated  by  a  circle.  Secondary 
modes  occur  where  an  unstable  beam  mode  (straight  line)  intersects  a  waveguide  mode 
(curved  line)  other  than  the  primary  mode.  These  are  indicated  by  squares.  Figure  2 
shows  the  temporal  growth  rate  of  the  primary  instability,  /m(u;/c)  plotted  versus  wave 
number  for  the  parameters  of  Fig.  1,  but  with  beam  current  It  =  2  kA.  These  figures, 
which  are  obtained  by  numerically  solving  Eq.  (23)  of  Ref.  5,  show  the  periodic  nature  of 
the  instability. 


2 


II.  Numerical  Results 


Simulations  were  performed  using  the  ELBA^  code.  ELBA  is  a  three-dimensional 
particle  code  which  simulates  a  beam  propagating  within  a  cylindrical  metallic  pipe.  The 
full  set  of  Maxwell’s  equations  along  with  the  full  relativistic  motion  of  the  beam  particles 
are  included.  The  ELBA  code  contains  numerous  beam  diagnostics,  including  an  rms 
emittance  diagnostic®  that  was  developed  for  beams  with  x-y  coupling. 

For  each  simulation,  initial  particle  positions  and  velocities  were  “matched”  to  the 
FODO  field  in  the  envelope  sense  as  in  Ref.  9,  including  an  ad  hoc  correction  for  space 
charge.  For  comparison  to  the  ELBA  runs,  simulations  were  also  performed  using  the 
SST^®  code.  The  SST  code  is  a  single-slice  version  of  the  ELBA  code  which  is  valid  in 
the  linndt  that  the  beam  radius  is  small  compared  to  the  length  scales  of  the  beam  and  the 
external  fields  in  the  axial  direction.  Because  the  derivatives  versus  the  axial  direction  are 
dropped  in  this  approximation,  the  SST  code  cannot  exhibit  the  instability.  This  allowed 
us  to  distinguish  physical  effects  (current  loss,  emittance  growth,  etc.)  that  may  be  caused 
by  the  instability  from  those  that  arise  from  imperfections  in  the  matching  scheme  (see, 
e.g..  Ref.  11). 

Growth  rates  were  measured  by  analyzing  the  TE\\  mode,  for  which  the  and  Er 
components  may  be  “projected  out”  from  the  electromagnetic  spectrum  in  a  straightfor¬ 
ward  manner.  The  growth  rate  is  then  obtained  as  a  function  of  u;  and  k  for  those  cases  in 
which  an  unstable  mode  grew  above  the  background  “noise” .  Because  our  simulation  takes 
place  in  a  coordinate  system  that  moves  with  the  beam  (r,0,(  =  ct  —  z),  these  growth 
rates  are  T/c  =  Im{u}/c  —  fc).  This  corresponds  to  the  theoretical  result  only  in  the  case 
that  Im{k)  =  0  as  was  assumed  in  numerically  solving  the  dispersion  relation. 

Simulations  were  performed  for  parameters  typical  of  a  high-current,  induction  accel¬ 
erator  beam,  /^  =  1  —  2  kA,  7o  =  7  and  normalized  rms  emittance,  en,rms  =  0-05  —  0.16 
cm-rad.  Simulation  parameters  are  listed  in  Table  I.  In  a  typical  run,  a  beam  of  length 
Z/fc  =  250  cm  was  transported  over  10  meters  in  the  presence  of  the  external  focusing  fields. 
It  was  not  practical  to  simulate  a  case  in  which  the  primary  mode  is  “cut-off,”  as  given  by 
Eq.  (2).  For  the  parameters  of  Fig.  2,  for  instance,  this  would  require  either  7  <  1.4  or 
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Tg  <  0.43  cm. 

In  each  run,  we  found  a  growing  TEu  mode  as  predicted  by  theory,  but  without 
a  growing  displacement  of  the  beam  centroid.  In  fact,  the  instability  saturated  at  a  low 
level,  apparently  due  to  mode  competition  between  the  various  wave  numbers  and,  possibly, 
between  the  primary  and  secondary  modes. 

As  an  example,  we  return  to  the  parameters  of  Fig.  2.  Here,  we  simulated  an  Lf,  — 
250  cm,  /fc  =  2  kA,  7  =  7,en,rTns  =  0.05  cm-rad  beam  propagating  for  10  meters.  The 
following  was  observed: 

1)  As  the  TEii  mode  grew,  the  peak  frequency  of  the  spectrum  shifted  continuously. 
Exponential  growth  rates  could  only  be  obtained  by  plotting  the  log  of  the  peak  of  the 
Fourier  transform  of  the  TEu  mode  versus  propagation  distance  z,  as  in  Fig.  3.  The 
points  in  Fig.  3  show  approximately  linear  growth  over  the  first  300  cm  of  propagation 
with  growth  rate  F/c  =  0.023  cm“^  and  frequency  tv/c  =  0.58  ±  0.12  cm“^.  Theoretical 
vsdues  for  this  case  are  F/c  =  0.024  cm“^  and  frequency  u>/c  =  0.67 ±0.04  cm~\  where  we 
are  quoting  the  half- width  at  half-maximum  of  the  F/c  versus  w/c  curve  as  the  expected 
frequency  spread. 

2)  Early  in  the  simulation  the  wave  number  spectrum,  shown  in  Fig.  4,  resembled 
that  expected  from  Fig.  2.  At  later  times  the  spectrum  is  broadband.  This  is  shown  in 
Fig.  5. 

3)  No  net  displacement  of  the  beam  centroid  was  observed.  In  fact,  the  instability 
had  no  discernable  effect  on  the  macroscopic  parameters  of  the  beam.  One  parameter  that 
did  evolve  over  the  course  of  the  simulation  was  the  emittance,  which  grew  as  a  result  of 
imperfections  in  the  initial  match  of  the  beam  to  the  external  fields.  This  growth  was  from 
Cn.rmj  =  0-05  cm-rad  to  0.20  cm-rad.  In  an  SST  run  with  identical  parameters,  emittance 
grew  from  en.rm*  =  0-05  cm-rad  to  0.13  cm-rad.  It  is  not  clear  whether  the  difference 
in  emittance  growth  between  the  ELBA  and  SST  runs  is  caused  by  the  instability  or  by 
differing  particle  statistics  between  the  two  codes.  The  latter  is  likely,  however,  as  the  final 
emittance  in  the  SST  code  was  reduced  somewhat  (<  10  %)  when  the  number  of  particles 
was  doubled. 

Results  for  all  runs  are  summarized  in  Table  II.  Growth  rates  and  frequencies  were 
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determined  as  described  above  for  the  Bqkq  =  1000  G/cm,  Vg  =  3.0  cm  case,  with  the¬ 
oretical  values  being  taken  from  numerical  solutions  of  Eq.  (23)  of  Ref.  5.  The  two 
runs  at  Bgkq  =  200  G/cm  showed  growth  rates  that  were  too  small  to  measure  reliably. 
The  Bqkq  =  1000  G/cm,  rg  =  4.0  cm  case,  run  D,  was  interesting  in  that  the  peak  fre¬ 
quency  of  the  TEii  mode  spectrum  clearly  jumped  between  the  primary  instability  at 
tjjjc  —  0.50  ±  0.03  cm~^,  and  a  higher  frequency  loJc  —  0.64  ±  0.02  cm“^.  The  value  listed 
in  Table  II  for  this  case,  w/c  =  0.54  ±  0.07  cm~^,  reflects  a  weighted  average  of  these 
two  frequencies.  The  higher  frequency  appears  to  correspond  to  the  secondary  mode  at 
u)/c  =  0.67  cm~^.  From  numerical  solutions  of  the  dispersion  relation,  the  secondary  mode 
should  have  a  growth  rate  T/c  <  0.002  cm~^  and,  therefore,  was  not  expected  to  play  a 
role  in  these  simulations. 
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III.  Conclusions 


An  electromagnetic  instability  on  an  intense  beam  in  a  FODO  lattice  was  studied 
via  numerical  simulation.  We  found  evidence  of  the  instability,  with  growth  rates  and 
frequencies  in  reasonable  agreement  with  theory.  The  inbi,...bility,  however,  was  found  to 
saturate  at  a  low  level,  with  field  amplitudes  two  orders  of  magnitude  lower  than  those 
observed  in  the  helical  quadrupole  case.  It  had  no  discernable  effect  on  the  macroscopic 
properties  of  the  beam. 

This  result  is  in  sharp  contrast  to  the  corresponding  instability  in  helical  quadrupole 
fields.  With  helical  quadrupole  fields,  the  instability  has  been  observed  in  simulations^’® 
and  in  experiments^  to  be  extremely  disruptive  to  beam  transport. 

We  also  find  that  the  instability  is  broad-band  in  both  frequency  and  wave  number.  We 
speculate  that  significant  mode  competition  is  occurring  both  among  the  various  unstable 
wave  numbers  of  the  primary  instability  and,  possibly,  between  the  primary  instability  and 
the  higher  frequency  secondary  instability.  It  appears  that  the  secondary  instability  may 
play  a  role  in  the  growth  of  the  primary  instability  despite  the  fact  that  its  theoretical 
growth  rate  is  more  than  an  order  of  magnitude  below  that  of  the  primary  instability. 
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Table  1 


Run 

h  (kA) 

kqBq  (G/cm) 

kq  (cm 

1)  Tg  (cm) 

A 

1.0 

200. 

0.60 

3.0 

B 

1.0 

200. 

0.40 

3.0 

C 

2.0 

1000. 

0.55 

3.0 

D 

2.0 

1000. 

0.55 

4.0 

Simulation  parameters.  Other  beam  parameters  '.vere  7  = 

7  and  en,rms  =0.16  cm-rad 

for  the  If,  = 

1  kA  cases  and  €n,rmt 

=  0.05  cm-rad  for  the  /j,  = 

2  kA  cases. 
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Table  II 


THEORY 


SIMULATION 


Run 

r/c 

a>/ c 

r/c 

ijjjc 

A 

0.0057 

0.62  ±  0.01 

- 

0.63  ±  0.03 

B 

0.0069 

0.71  ±0.02 

- 

0.55  ±0.07 

C 

0.024 

0.67  ±  0.03 

0.023 

0.58  ±0.12 

D 

0.021 

0.50  ±  0.01 

0.027 

0.54  ±  0.07 

Theoretical  values  of  peak  linear  growth  rate  T/c  and  frequency  u>/c  and  corresponding 
simulation  results.  Simulation  growth  rates  for  runs  A  and  B  could  not  be  determined 
reliably  because  the  TE\\  mode  was  largely  masked  by  particle  “noise”  in  those  cases. 
Growth  rates  and  frequencies  are  in  units  of  cm“^. 
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Fig.  2.  TemporaJ  growth  rate  of  the  primary  mode  for  <*;  >  0  as  a  function  of  wave  number  k 
for  parameters  listed  in  Table  I,  Case  C. 
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12 


Fig.  3.  The  log  of  the  peak  of  the  Fourier  transform,  in  frequency  space,  of  the  component 
of  the  TEn  mode  versus  propagation  distance  for  the  Bgkg  =  1000  G/cm,  Vg  =  3.0  cm 
case  (Case  C  in  Tables  I  and  11).  A  linear  least-squares  fit  to  the  first  nine  points  is 
also  plotted. 
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Fig.  4.  Wave  number  spect 
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